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Summary. The theory of cosmological perturbations has become a cornerstone of 
modern quantitative cosmology since it is the framework which provides the link 
between the models of the very early Universe such as the inflationary Universe 
scenario (which yield causal mechanisms for the generation of fluctuations) and 
the wealth of recent high-precision data on the spectrum of density fluctuations 
and cosmic microwave anisotropies. In these lectures, I provide an overview of the 
classical and quantum theory of cosmological fluctuations. 

Crucial points in both the current inflationary paradigm ^ [2] of the early 
Universe and in proposed alternatives such as the Pre-Big-Bang [3] and Ekpyrotic 
|3] scenarios are that, first, the perturbations are generated on microscopic scales 
as quantum vacuum fluctuations, and, second, that via an accelerated expansion of 
the background geometry (or by a contraction of the background), the wavelengths 
of the fluctuations become much larger than the Hubble radius for a long period 
of cosmic evolution. Hence, both Quantum Mechanics and General Relativity are 
required in order to understand the generation and evolution of fluctuations. 

As a guide to develop the physical intuition for the evolution of inhomogeneities, 
I begin with a discussion of the Newtonian theory of fluctuations, applicable at 
late times and on scales smaller than the Hubble radius. The analysis of super- 
Hubble fluctuations requires a general relativistic analysis. I first review the classical 
relativistic theory of fluctuations, and then discuss their quantization. I conclude 
with a brief overview of two applications of the theory of cosmological fluctuations: 
the trans-Planckian "problem" of inflationary cosmology and the current status 
of the study of the back-reaction of cosmological fluctuations on the background 
space-time geometry. Most of this article is based on the review to which the 
reader is referred to for the details omitted in these lecture notes. 



1 Motivation 

As described in the lectures by Tegmark at this school observational 
cosmology is currently in its golden years. Using a variety of observational 
techniques, physicists and astronomers are exploring the large-scale structure 
of the Universe. The Cosmic Microwave Background (CMB) is the observa- 
tional window which in recent years has yielded the most information. The 
anisotropies in the CMB have now been detected on a wide range of angu- 
lar scales, giving us a picture of the Universe at the time of recombination, 
the time that the cosmic photons last scattered. Large-scale galaxy redshift 
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surveys are providing us with increasingly accurate power spectra of the dis- 
tribution of objects in the Universe which emit light, which - modulo the 
question whether light in fact traces mass (this is the issue of the cosmic 
bias) - gives us the distribution of mass at the present time. Analyses of 
the spectra of quasar absorption line systems and weak gravitational lensing 
surveys are beginning to give us complementary information about the dis- 
tribution of matter (independent of whether this matter in fact emits light, 
thus shedding light on the biasing issue). The analysis of weak gravitational 
lensing maps is in fact sensitive not only to the baryonic but also to the dark 
matter, and promises to give a technique which unambiguously reveals where 
the dark matter is concentrated. X-ray telescopes are providing additional 
information on the distribution of sources which emit X-rays. 

The current data fits astonishingly well with the current paradigm of early 
Universe cosmology, the inflationary Universe scenario pp. However, it is im- 
portant to keep in mind that what is tested observationally is the paradigm 
that the primordial spectrum of inhomogeneities was scale-invariant and pre- 
dominantly adiabatic (these terms will be explained in the following Section) , 
and that there might exist other scenarios of the very early Universe which 
do not yield inflation but predict a scale-invariant adiabatic spectrum. For 
example, within both the Pre-Big-Bang |2j and the Ekpyrotic scenarios 0] 
there may be models which yield such a spectrum 1 One should also not for- 
get that topological defect models of structure formation (see e.g. [r31l2T)l I2T] 
for reviews) naturally yield a scale-invariant spectrum, however of primor- 
dial isocurvature nature and thus no longer compatible with the latest CMB 
anisotropy results. 

The theory of cosmological perturbations is what allows us to connect 
theories of the very early Universe with the data on the large-scale structure 
of the Universe at late times and is thus of central importance in modern cos- 
mology. The techniques discussed below are applicable to most scenarios of 
the very early Universe. Most specific applications mentioned, however, will 
be within the context of the inflationary Universe scenario. To understand 
what the key requirements for a viable theory of cosmological perturbations 
are, recall the basic space-time diagram for inflationary cosmology (Figure 
1): Since, during the phase of standard cosmology £_r < t < t , where tn cor- 

1 Note, however, that whereas the simplest inflationary models yield an almost 
scale-invariant n — 1 spectrum of fluctuations, as discussed in detail in these 
lectures, this is not the case for the simplest models of Pre-Big-Bang type nor 
for four dimensional descriptions of the Ekpyrotic scenario. In the case of single 
field realizations of Pre-Big-Bang cosmology, a spectrum with spectral index 
n = 4 emerges [7 . In Ekpyrotic cosmology, the value of the index of the final 
power spectrum is under active debate. Most studies conclude either that the 
spectral index is n = 3 |81l§l ll(Jlllllll2| . or that the result is ill-defined because of 
the singularities at the bounce 1 131 II 4| (see, however, (151 1161 HT) for arguments 
in support of a final scale- invariant spectrum) . See also |18| for criticisms of the 
basic setup of the Ekpyrotic scenario. 
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Fig. 1. Space-time diagram (sketch) showing the evolution of scales in inflationary 
cosmology. The vertical axis is time, and the period of inflation lasts between U 
and tn, and is followed by the radiation-dominated phase of standard big bang 
cosmology. During exponential inflation, the Hubble radius H^ 1 is constant in 
physical spatial coordinates (the horizontal axis), whereas it increases linearly in 
time after tn. The physical length corresponding to a fixed comoving length scale 
labelled by its wavenumber k increases exponentially during inflation but increases 
less fast than the Hubble radius (namely as t 1 ^ 2 ), after inflation. 

responds to the end of inflation, and to denotes the present time, the Hubble 
radius Zjy(t) = i7 _1 (t) expands faster that the physical wavelength associ- 
ated with a fixed comoving scale, the wavelength becomes larger than the 
Hubble radius as we go backwards in time. However, during the phase of ac- 
celerated expansion (inflation), the physical wavelength increases much faster 
than the Hubble radius, and thus at early times the fluctuations emerged at 
micro-physical sub-Hubble scales. The idea is that micro-physical processes 
(as we shall see, quantum vacuum fluctuations) are responsible for the origin 
of the fluctuations. However, during the period when the wavelength is super- 
Hubblc, it is essential to describe the fluctuations using General Relativity. 
Thus, both Quantum Mechanics and General Relativity are required to suc- 
cessfully describe the generation and evolution of cosmological fluctuations. 

A similar conclusion can be reached when considering the space-time di- 
agram in a model of Pre-Big-Bang or Ekpyrotic type, where the Universe 
starts out in a contracting phase during which the Hubble radius contracts 
faster than the physical length corresponding to a fixed comoving scale (see 
Figure 2) . The contracting phase ends at a cosmological bounce, after which 
the Universe is assumed to follow the same evolution history as it does in 
standard Big Bang cosmology. As in inflationary cosmology, quantum vac- 
uum fluctuation on sub-Hubble scales (in this case in the contracting phase) 
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Fig. 2. Space-time diagram (sketch) showing the evolution of scales in a cosmology 
of PBB or Ekpyrotic type. The axes are as in Figure 1. Times earlier than ts 
correspond to the contracting phase, times after describe the post-bounce phase of 
expansion as described in standard cosmology. The Hubble radius decreases relative 
to a fixed comoving scale during the contracting phase, and increases faster in the 
expanding phase. Fluctuations of cosmological interest today are generated sub- 
Hubble but propagate super-Hubble for a long time interval. 

are assumed to be the seeds of the inhomogeneities observed today. For a 
long time period, the scale of the fluctuation is super-Hubble. 

Thus, we see that in inflationary cosmology as well as in Pre-Big-Bang 
and Ekpyrotic-type models, both Quantum Mechanics and General Relativ- 
ity are required to understand the generation and evolution of cosmological 
perturbations. 

2 Newtonian Theory of Cosmological Perturbations 
2.1 Introduction 

The growth of density fluctuations is a consequence of the purely attractive 
nature of the gravitational force, fmaginc (first in a non-expanding back- 
ground) a density excess 5p localized about some point x in space. This 
fluctuation produces an attractive force which pulls the surrounding mat- 
ter towards x. The magnitude of this force is proportional to Sp. Hence, by 
Newton's second law 

Sp ~ GSp, (1) 

where G is Newton's gravitational constant. Hence, there is an exponential 
instability of flat space-time to the development of fluctuations. 
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Obviously, in General Relativity it is inconsistent to consider density fluc- 
tuations in a non-expanding background. If we consider density fluctuations 
in an expanding background, then the expansion of space leads to a friction 
term in (Q. Hence, instead of an exponential instability to the development 
of fluctuations, the growth rate of fluctuations in an expanding Universe will 
be as a power of time. It is crucial to determine what this power is and how 
it depends both on the background cosmological expansion rate and on the 
length scale of the fluctuations. 

We will be taking the background space-time to be homogeneous and 
isotropic, with a metric given by 

ds 2 = dt 2 - a(t) 2 dx 2 , (2) 

where t is physical time, dx 2 is the Euclidean metric of the spatial hyper- 
surfaces (here taken for simplicity to be spatially flat), and a(t) denoting the 
scale factor, in terms of which the expansion rate is given by H(t) = a/a. The 
coordinates x used above are "comoving" coordinates, coordinates painted 
onto the expanding spatial hypersurfaces. Note, however, that in the following 
two subsections x will denote the physical coordinates, and q the comoving 
ones. 

The materials covered in this section are discussed in several excellent 
textbooks on cosmology, e.g. in |221 1231 I25| . 



2.2 Perturbations about Minkowski Space-Time 

To develop some physical intuition, we first consider the evolution of hydro- 
dynamical matter fluctuations in a fixed non-expanding background. Note 
that in this case the background Einstein equations are not satisfied. 

In this context, matter is described by a perfect fluid, and gravity by 
the Newtonian gravitational potential ip. The fluid variables are the energy 
density p, the pressure p, the fluid velocity v, and the entropy density S. The 
basic hydrodynamical equations are 

p + V • (pv) = 

v + (v • V)v + - X7p + X7tp = 
P 

VV = ^Gp (3) 
S+ (v V)5 = 

P = p(p, s) . 

The first equation is the continuity equation, the second is the Euler (force) 
equation, the third is the Poisson equation of Newtonian gravity, the fourth 
expresses entropy conservation, and the last describes the equation of state 
of matter. The derivative with respect to time is denoted by an over-dot. 



6 Robert H. Brandenberger 



The background is given by the background energy density p Q1 the back- 
ground pressure po, vanishing velocity, constant gravitational potential ipo 
and constant entropy density Sq. As mentioned above, it does not satisfy the 
background Poisson equation. 

The equations for cosmological perturbations are obtained by perturbing 
the fluid variables about the background, 

P = Po + dp 
v = <5v 

P = Po + Sp (4) 

(fi = IPQ+Slfi 

S = S + 6S, 

where the fluctuating fields Sp, <5v, Sp, Sip and SS are functions of space and 
time, by inserting these expressions into the basic hydrodynamical equations 
(|2J), by linearizing, and by combining the resulting equations which are of 
first order in time to obtain the following second order differential equations 
for the energy density fluctuation bp and the entropy perturbation SS 

Sp - c 2 s V 2 5p - AirGpoSp = aV 2 5S (5) 
<55 =0, 

where the variables c 2 and a describe the equation of state 

8p = c 2 Sp + aSS (6) 

with 

4 - 0u M 

denoting the square of the speed of sound. 

What can we learn from these equations? First of all, since the equations 
are linear, we can work in Fourier space. Each Fourier component 5pk(t) of 
the fluctuation field <5p(x, t) 

S P (x,t) = J e^'Spkit) (8) 

evolves independently. 

There are various types of fluctuations. If the entropy fluctuation 5S van- 
ishes, we have adiabatic fluctuations. If the entropy fluctuation SS is non- 
vanishing but dp = 0, we speak on an entropy fluctuation. 

The first conclusions we can draw from the basic perturbation equations 
(0 are that 

1) entropy fluctuations do not grow, 

2) adiabatic fluctuations are time-dependent, and 

3) entropy fluctuations seed an adiabatic mode. 
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Taking a closer look at the equation of motion (J5J) for Sp, we see that the 
third term on the left hand side represents the force due to gravity, a purely 
attractive force yielding an instability of flat space-time to the development 
of density fluctuations (as discussed earlier, see Q). The second term on the 
left hand side of JSJ) represents a force due to the fluid pressure which tends 
to set up pressure waves. In the absence of entropy fluctuations, the evolution 
of Sp is governed by the combined action of both pressure and gravitational 
forces. 

Restricting our attention to adiabatic fluctuations, we see from (|3J that 
there is a critical wavelength, the Jeans length, whose wavenumber kj is given 

by 

kj=(^) 1/2 . (9) 

Fluctuations with wavelength longer than the Jeans length (k <C kj) grow 
exponentially 

8p k {t) ~ e" kt with Lo k ~ 4{irGpo) 1/2 (10) 

whereas short wavelength modes (fc 3> kj) oscillate with frequency Uk ~ c s k. 
Note that the value of the Jeans length depends on the equation of state of 
the background. For a background dominated by relativistic radiation, the 
Jeans length is large (of the order of the Hubble radius _ff _1 (i)), whereas for 
pressure-less matter the Jeans length goes to zero. 



2.3 Perturbations about an Expanding Background 

Let us now improve on the previous analysis and study Newtonian cosmologi- 
cal fluctuations about an expanding background. In this case, the background 
equations are consistent (the non-vanishing average energy density leads to 
cosmological expansion). However, we are still neglecting general relativistic 
effects (the fluctuations of the metric). Such effects turn out to be dominant 
on length scales larger than the Hubble radius H~ l (t), and thus the analysis 
of this section is applicable only to scales smaller than the Hubble radius. 

The background cosmological model is given by the energy density po(t), 
the pressure po(t), and the recessional velocity vo = H(t)x, where x is the 
Euclidean spatial coordinate vector ( "physical coordinates" ) . The space- and 
timc-dcpcndent fluctuating fields are defined as in the previous section: 

p(t,x) =po(t)(l + S e (t,x)) 

v(t,x) = v (t,x) + <Jv(t,x) (11) 
p(t,x) = p (t) + 5p(t,x) , 

where S c is the fractional energy density perturbation (we are interested in 
the fractional rather than in the absolute energy density fluctuation!), and 
the pressure perturbation Sp is defined as in ©. In addition, there is the 
possibility of a non- vanishing entropy perturbation defined as in jl)!. 
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We now insert this ansatz into the basic hydrodynamical equations (J3J, 
linearize in the perturbation variables, and combine the first order differential 
equations for <5 e and Sp into a single second order differential equation for Sp e . 
The result simplifies if we work in "comoving coordinates" q which are the 
coordinates painted onto the expanding background, i.e. 

x(t) = o(t)q(t) . (12) 

After a substantial amount of algebra, we obtain the following equation which 
describes the time evolution of density fluctuations: 

<? a 
S e + 2H6 e - ^% ~ InGpoSe = ^6S , (13) 

where the subscript q on the V operator indicates that derivatives with re- 
spect to comoving coordinates are used. In addition, we have the equation of 
entropy conservation 

6'S = . (14) 

Comparing with the equations (jSJ obtained in the absence of an expand- 
ing background, we see that the only difference is the presence of a Hubble 
damping term in the equation for 5 e . This term will moderate the exponen- 
tial instability of the background to long wavelength density fluctuations. In 
addition, it will lead to a damping of the oscillating solutions on short wave- 
lengths. More specifically, for physical wavenumbers k p <C kj (where kj is 
again given by ©), and in a matter-dominated background cosmology, the 
general solution of in the absence of any entropy fluctuations is given by 

S k (t) = cit 2 ' 3 + cat -1 , (15) 

where c\ and ci are two constants determined by the initial conditions, and 
we have dropped the subscript e in expressions involving S e . There are two 
fundamental solutions, the first is a growing mode with Sk(t) ~ a(t), the 
second a decaying mode with Sk(t) ~ t . On short wavelength, one obtains 
damped oscillatory motion: 

S k {t) ~ a - 1/2 (t)exp(±ic s k J dt'a" 1 ^)) . (16) 

As a simple application of the Newtonian equations for cosmological per- 
turbations derived above, let us compare the predicted cosmic microwave 
background (CMB) anisotropics in a spatially flat Universe with only bary- 
onic matter - Model A - to the corresponding anisotropics in a flat Uni- 
verse with mostly cold dark matter (pressure-less non-baryonic dark matter) 
- Model B. We start with the observationally known amplitude of the relative 
density fluctuations today (time to), and we use the fact that the amplitude of 
the CMB anisotropies on the angular scale 9{k) corresponding to the comov- 
ing wavenumber k is set by the value of the primordial gravitational potential 
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<p - introduced in the following section - which in turn is related to the value 
of the primordial density fluctuations at Hubble radius crossing (and not to 
its value of the time t rec ). See e.g. Chapter 17 of 

In Model A, the dominant component of the pressure-less matter is cou- 
pled to radiation between t eq and t rec , the time of last scattering. Thus, the 
Jeans length is comparable to the Hubble radius. Therefore, for comoving 
galactic scales, k 3> kj in this time interval, and thus the fractional density 
contrast decreases as a{t)^ x ^ 2 . In contrast, in Model B, the dominant com- 
ponent of pressure-less matter couples only weakly to radiation, and hence 
the Jeans length is negligibly small. Thus, in Model B, the relative density 
contrast grows as a(t) between t eq and t rec . In the time interval t rec < t < to, 
the fluctuations scale identically in Models A and B. Summarizing, we con- 
clude, working backwards in time from a fixed amplitude of 8k today, that the 
amplitudes of 8k(t eq ) in Models A and B (and thus their primordial values) 
are related by 



Hence, in Model A (without non-baryonic dark matter) the CMB anisotropies 
are predicted to be a factor of about 30 larger jJH] than in Model B, way 
in excess of the recent observational results. This is one of the strongest 
arguments for the existence of non-baryonic dark matter. 

2.4 Characterizing Perturbations 

Let us consider perturbations on a fixed comoving length scale given by a co- 
moving wavenumber k. The corresponding physical length increases as a(t). 
This is to be compared to the Hubble radius i? _1 (t) which scales as t pro- 
vided ait) grows as a power of t. In the late time Universe, a(t) ~ t 1 / 2 in the 
radiation-dominated phase (i.e. for t < t eq , and a(t) ~ t 2 ^ 3 in the matter- 
dominated period (t eq < t < to). Thus, we see that at sufficiently early times, 
all comoving scales had a physical length larger than the Hubble radius. If 
we consider large cosmological scales (e.g. those corresponding to the ob- 
served CMB anisotropies or to galaxy clusters), the time tn(k) of "Hubble 
radius crossing" (when the physical length was equal to the Hubble radius) 
was in fact later than t eq . As we will see in later sections, the time of Hub- 
ble radius crossing plays an important role in the evolution of cosmological 
perturbations. 

Cosmological fluctuations can be described either in position space or 
in momentum space. In position space, we compute the root mean square 
mass fluctuation 8M/M(k,t) in a sphere of radius I — 2n/k at time t. A 
scale-invariant spectrum of fluctuations is defined by the relation 




(17) 



— (k,t H (k)) = const. 



(18) 
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Such a spectrum was first suggested by Harrison [27J and Zeldovich [2H| as 
a reasonable choice for the spectrum of cosmological fluctuations. We can 
introduce the "spectral index" n of cosmological fluctuations by the relation 

( 6 -M-) 2 (k,t H (k)) ~ k»-\ (19) 

and thus a scale-invariant spectrum corresponds to n = 1. 

To make the transition to the (more frequently used) momentum space 
representation, we Fourier decompose the fractional spatial density contrast 

zk-x 



5 e (x,i) = / d d fc<5 e (k,i)e lkx . (20) 
The power spectrum Pg of density fluctuations is defined by 

P s (k) = k 3 \6~ e (k)\ 2 , (21) 

where k is the magnitude of k, and we have assumed for simplicity a Gaussian 
distribution of fluctuations in which the amplitude of the fluctuations only 
depends on k. 

We can also define the power spectrum of the gravitational potential ip: 

P v (k) = k 3 \6~<p(k)\ 2 . (22) 

These two power spectra are related by the Poisson equation iffil 

P v (k) ~ k~ 4 P s (k). (23) 

In general, the condition of scale-invariance is expressed in momentum 
space in terms of the power spectrum evaluated at a fixed time. To obtain 
this condition, we first use the time dependence of the fractional density 
fluctuation from l|15l) to determine the mass fluctuations at a fixed time 
t > tu{k) > t eq (the last inequality is a condition on the scales considered) 

^ M = ^ 1 (24) 
The time of Hubble radius crossing is given by 

a{t H {k))k- x = 2t H {k), (25) 

and thus 

t H {kf/ 2 ~ k- 1 . (26) 
Inserting this result into l|24|l making use of (|19Jl we find 



(27) 
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Since, for reasonable values of the index of the power spectrum, SM/M(k,t) 
is dominated by the Fourier modes with wavenumber k, we find that l|27[l 
implies 

|£| 2 ~fc", (28) 

or, equivalently, 

P v (k) ~ k n - 1 . (29) 
2.5 Matter Fluctuations in the Radiation Era 

Let us now briefly consider fluctuations in the radiation dominated epoch. 
We are interested in both the fluctuations in radiation and in matter (cold 
dark matter). In the Newtonian treatment, Eq. (|13|1 is replaced by separate 
equations for each matter fluid component (these components are designated 
by the labels A or B): 

6 A + 2Kb A - v 2 A a- 2 V 2 5 A = AttG^PbSb , (30) 

B 

where ps indicate the background densities, and 5b the fractional density 
fluctuations. The velocities of the respective fluid components are denoted 
by Vb, with v 2 . = 1/3 for radiation and v m — for cold dark matter. 

In the radiation dominated epoch, the evolution of the fluctuations in 
radiation is to a first approximation (in the ratio of the background densities) 
independent of the cold matter content. Inserting the expansion rate for this 
epoch, we thus immediately obtain 

5 r (t) ~ a{tf (31) 

on scales much larger than the Hubble scale, i.e. k <C kjj , whereas S r under- 
goes damped oscillatory motion on smaller scales. 

The evolution of the matter fluctuation 5 m is more complicated. Its equa- 
tion of motion is dominated by the source term coming from S r . What results 
is logarithmic growth of the amplitude of S m , instead of the growth propor- 
tional to a(t) which would occur on these scales in the absence of radiation. 
This damping effect on matter fluctuations due to the presence of radiation 
is called the "Meszaros effect". It leads to a turnover in the spectrum of 
cosmological fluctuations at a scale k eq which crosses the Hubble radius at 
the time of equal matter and radiation. On larger scales (k < k eq ), one has 
the primordial power spectrum with spectral index n, on smaller scales, to a 
first approximation, the spectral index changes to n — 4. The details of the 
power spectrum on small scales depend largely on the specifics of the matter 
content in the Universe. One can write 



Pfinal(k,t) = T(k,t)P (k,t) 



(32) 
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where Pq is the primordial power spectrum extrapolated to late times with 
unchanged spectral index, and Pfmai denotes the actual power spectrum 
which depends on effects such as the ones mentioned above. For more details 
see e.g. 



3 Relativistic Theory of Cosmological Fluctuations 
3.1 Introduction 

The Newtonian theory of cosmological fluctuations discussed in the previ- 
ous section breaks down on scales larger than the Hubble radius because it 
neglects perturbations of the metric, and because on large scales the metric 
fluctuations dominate the dynamics. 

Let us begin with a heuristic argument to show why metric fluctuations 
are important on scales larger than the Hubble radius. For such inhomo- 
geneities, one should be able to approximately describe the evolution of 
the space-time by applying the first Friedmann-Lemaitre-Robertson- Walker 
(FLRW) equation of homogeneous and isotropic cosmology to the local Uni- 
verse (this approximation is made more rigorous in |29|): 

,au2 8itG 

(-) =— P- 03) 

Based on this equation, a large-scale fluctuation of the energy density will 
lead to a fluctuation ("5a") of the scale factor a which grows in time. This is 
due to the fact that self gravity amplifies fluctuations even on length scales 
A greater than the Hubble radius. 

This argument is made rigorous in the following analysis of cosmological 
fluctuations in the context of general relativity, where both metric and matter 
inhomogeneities are taken into account. We will consider fluctuations about 
a homogeneous and isotropic background cosmology, given by the metric , 
which can be written in conformal time r\ (defined by dt — a(t)drj) as 

ds 2 = a(r]) 2 (dr] 2 ~ dx 2 ) . (34) 

The evolution of the scale factor is determined by the two FLRW equations, 
and 

p=-3H(p + p), (35) 

which determine the expansion rate and its time derivative in terms of the 
equation of state of the matter, whose background stress-energy tensor can 
be written as 

( p \ 
-p 
-p 
\0 -pj 



TZ = 



(36) 
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The theory of cosmological perturbations is based on expanding the Ein- 
stein equations to linear order about the background metric. The theory was 
initially developed in pioneering works by Lifshitz 30 . Significant progress in 
the understanding of the physics of cosmological fluctuations was achieved by 
Bardeen [31] who realized the importance of subtracting gauge artifacts (see 
below) from the analysis (see also [22] )• The following discussion is based on 
Part I of the comprehensive review article Other reviews - in some cases 
emphasizing different approaches - are [23 Ell OH CM! • 

3.2 Classifying Fluctuations 

The first step in the analysis of metric fluctuations is to classify them ac- 
cording to their transformation properties under spatial rotations. There are 
scalar, vector and second rank tensor fluctuations. In linear theory, there is 
no coupling between the different fluctuation modes, and hence they evolve 
independently (for some subtleties in this classification, see EZ]). 



We begin by expanding the metric about the FLRW background metric 
g$ given by 



The background metric depends only on time, whereas the metric fluctuations 
Sg^v depend on both space and time. Since the metric is a symmetric tensor, 
there are at first sight 10 fluctuating degrees of freedom in Sg^. 

There are four degrees of freedom which correspond to scalar metric fluc- 
tuations (the only four ways of constructing a metric from scalar functions): 



where the four fluctuating degrees of freedom are denoted (following the no- 
tation of U2|) 0, -B, E, and ip, a comma denotes the ordinary partial derivative 
(if we had included spatial curvature of the background metric, it would have 
been the covariant derivative with respect to the spatial metric), and Sij is 
the Kronecker symbol. 

There are also four vector degrees of freedom of metric fluctuations, con- 
sisting of the four ways of constructing metric fluctuations from three vectors: 



where Si and Fi are two divergence-less vectors (for a vector with non- 
vanishing divergence, the divergence contributes to the scalar gravitational 
fluctuation modes). 

Finally, there are two tensor modes which correspond to the two polar- 
ization states of gravitational waves: 



9v.i> = g, 




(37) 




(38) 
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^ = -° 2 (o4)' (40) 

where hij is trace-free and divergence-less 

K = K{j=U- (41) 

Gravitational waves do not couple at linear order to the matter fluctua- 
tions. Vector fluctuations decay in an expanding background cosmology and 
hence are not usually cosmologically important. The most important fluctu- 
ations, at least in inflationary cosmology are the scalar metric fluctuations, 
the fluctuations which couple to matter inhomogeneities and which are the 
relativistic generalization of the Newtonian perturbations considered in the 
previous section. 

3.3 Gauge Transformation 

The theory of cosmological perturbations is at first sight complicated by the 
issue of gauge invariance (at the final stage, however, we will see that we can 
make use of the gauge freedom to substantially simplify the theory) . The co- 
ordinates t, x of space-time carry no independent physical meaning. They are 
just labels to designate points in the space-time manifold. By performing a 
small-amplitude transformation of the space-time coordinates (called "gauge 
transformation" in the following), we can easily introduce "fictitious" fluc- 
tuations in a homogeneous and isotropic Universe. These modes are "gauge 
artifacts" . 

We will in the following take an "active" view of gauge transforma- 
tion. Let us consider two space-time manifolds, one of them a homogeneous 
and isotropic Universe Mo, the other a physical Universe M with inhomo- 
geneities. A choice of coordinates can be considered to be a mapping V be- 
tween the manifolds Mo and M. Let us consider a second mapping T> which 
will map the same point (e.g. the origin of a fixed coordinate system) in Mo 
into different points in M. Using the inverse of these maps V and V, we can 
assign two different sets of coordinates to points in M . 

Consider now a physical quantity Q (e.g. the Ricci scalar) on M, and the 
corresponding physical quantity on Mo Then, in the first coordinate 
system given by the mapping T>, the perturbation SQ of Q at the point 
p G M is defined by 

SQ(p) = Q(p)-Q {0) (p- 1 (p)). (42) 

Analogously, in the second coordinate system given by T>, the perturbation 
is defined by 

S~Q(p) = Q{p)-Q (0) ip-\p)). (43) 

The difference 
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AQ(p) = 5~Q(p) - SQ( P ) (44) 

is obviously a gauge artifact and carries no physical significance. 

Some of the metric perturbation degrees of freedom introduced in the first 
subsection are gauge artifacts. To isolate these, we must study how coordinate 
transformations act on the metric. There are four independent gauge degrees 
of freedom corresponding to the coordinate transformation 

x» -> x>* ^x fl + ^. (45) 

The zero (time) component £° of £^ leads to a scalar metric fluctuation. The 
spatial three vector can be decomposed 

f = 4+7%- (46) 

(where 7 lJ is the spatial background metric) into a transverse piece £| r which 
has two degrees of freedom which yield vector perturbations, and the second 
term (given by the gradient of a scalar £) which gives a scalar fluctuation. 
To summarize this paragraph, there are two scalar gauge modes given by 
£° and £, and two vector modes given by the transverse three vector Q. r . 
Thus, there remain two physical scalar and two vector fluctuation modes. 
The gravitational waves are gauge-invariant. 

Let us now focus on how the scalar gauge transformations (i.e. the trans- 
formations given by £° and £) act on the scalar metric fluctuation variables 
4>, B, E, and ip. An immediate calculation yields: 

4> = 0--£°-(£ )' 

a 

B = B + £° - £' (47) 
E = E-i 

i> = ^ + 

a 

where a prime indicates the derivative with respect to conformal time rj. 

There are two approaches to deal with the gauge ambiguities. The first 
is to fix a gauge, i.e. to pick conditions on the coordinates which completely 
eliminate the gauge freedom, the second is to work with a basis of gauge- 
invariant variables. 

If one wants to adopt the gauge-fixed approach, there are many different 
gauge choices. Note that the often used synchronous gauge determined by 
Sg ^ = does not totally fix the gauge. A convenient system which com- 
pletely fixes the coordinates is the so-called longitudinal or conformal 
Newtonian gauge defined by B = E = 0. 

If one prefers a gauge-invariant approach, there are many choices of gauge- 
invariant variables. A convenient basis first introduced by [2] is the basis <P, \P 
given by 
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<P = <j> + -\(B - E')a\ (48) 
a 

& = tp- —(B — E') . (49) 
a 

It is obvious from the above equations that the gauge-invariant variables 
and coincide with the corresponding diagonal metric perturbations (j) an( i 
tp in longitudinal gauge. 

Note that the variables defined above are gauge-invariant only under lin- 
ear space-time coordinate transformations. Beyond linear order, the structure 
of perturbation theory becomes much more involved. In fact, one can show 
|38| that the only fluctuation variables which are invariant under all coordi- 
nate transformations are perturbations of variables which are constant in the 
background space-time. 



3.4 Equation of Motion 

We begin with the Einstein equations 

G M „ = 8nGT^ , (50) 

where G^ v is the Einstein tensor associated with the space-time metric g^„, 
and T^ v is the energy-momentum tensor of matter, insert the ansatz for 
metric and matter perturbed about a FLRW background (g^(?7), ip^°\rj)): 

c^„(x,77) = g$(r)) +Sg fiv (x,r)) (51) 
(p(x,r)) = <p (r)) +S<p(x,r}) , (52) 

(where we have for simplicity replaced general matter by a scalar matter field 
if) and expand to linear order in the fluctuating fields, obtaining the following 
equations: 

5G^ = 8ttG5T^ . (53) 

In the above, 5g^ v is the perturbation in the metric and Sip is the fluctuation 
of the matter field tp. 

Note that the components 8G^ and ST^ are not gauge invariant. If we 
want to use the gauge-invariant approach, we note [5] that it is possible to 
construct a gauge-invariant tensor SGi 9 ^ M via 

<5G^ )0 = SG^ + i^G'^iB-E') 

6G (gi)o _ SG o + ((0 ) G o _ i(o) G fc)( B _ E'),i (54) 

SG ( f t)l = 8G) + ( {0) G'/){B - E') , 

where ^0)G^ denote the background values of the Einstein tensor. Analo- 
gously, a gauge-invariant linearized stress-energy tensor ST^ 9 ^ M can be de- 
fined. In terms of these tensors, the gauge-invariant form of the equations of 
motion for linear fluctuations reads 
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6G$> = teGSTjg) . (55) 

If we insert into this equation the ansatz for the general metric and matter 
fluctuations (which depend on the gauge) , only gauge- invariant combinations 
of the fluctuation variables will appear. 

In a gauge-fixed approach, one can start with the metric in longitudinal 
gauge 

ds 2 = a 2 [(1 + 2cj))dT] 2 - (1 - 2^)j lJ dx i dx j ] (56) 

and insert this ansatz into the general perturbation equations i|53|) ■ The short- 
cut of inserting a restricted ansatz for the metric into the action and deriving 
the full set of variational equations is justified in this case. 

Both approaches yield the following set of equations of motion: 

- 3H(H$ + + V 2 tf = 47rGa 2 (5T (9l) 

(H$ + . = ^Ga 2 bT\ 9l) (57) 
[(2H +H 2 )$ + H$ + W" + 2HW'] 5) 

+\v 2 D5) - p ik D, kj = -^Ga 2 ST^\ 

where D = <P — W and 7i = a' /a. If we work in longitudinal gauge, then 
ST (gi) i = §T j $ = ^ and ^ = ^ 

The first conclusion we can draw is that if no anisotropic stress is present 
in the matter at linear order in fluctuating fields, i.e. 8TJ- = for i ^ j, then 
the two metric fluctuation variables coincide: 

<P = <P. (58) 

This will be the case in most simple cosmological models, e.g. in theories with 
matter described by a set of scalar fields with canonical form of the action, 
and in the case of a perfect fluid with no anisotropic stress. 

Let us now restrict our attention to the case of matter described in terms 
of a single scalar field ip with action 

S = j dW^H^'Va - Vfr)] (59) 

(where g denotes the determinant of the metric) and we expand the matter 
field as 

(p(x,i]) = <p (r]) +6cp(x,r)) (60) 

in terms of background matter ipo and matter fluctuation 5ip(x., rj), then in 
longitudinal gauge (|57|l reduce to the following set of equations of motion 
(making use of (|58[0 

V 2 <j> - m<p - (H + 2U 2 )<j) = 4nG(tp' 5(p + v'a 2 S(p) 

4 +H<j) = 4nGip' 6ip (61) 
4' + 3H<f> + (H + 2H 2 )<t> = 4irG(tp 6tp - v'a 2 S(p) , 
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where V denotes the derivative of V with respect to (p. These equations can 
be combined to give the following second order differential equation for the 
relativistic potential <j>: 

4' + 2 ( H - ) 4 ~ V 2 + 2 ( U' - U^r |0 = O. (62) 

Let us now discuss this final result for the classical evolution of cosmolog- 
ical fluctuations. First of all, we note the similarities with the equation l|13fl 
obtained in the Newtonian theory. The final term in (|62fl is the force due to 
gravity leading to the instability, the second to last term is the pressure force 
leading to oscillations (relativistic since we are considering matter to be a 
relativistic field), and the second term is the Hubble friction term. For each 
wavenumber there are two fundamental solutions. On small scales (fc > H), 
the solutions correspond to damped oscillations, on large scales (k < H) the 
oscillations freeze out and the dynamics is governed by the gravitational force 
competing with the Hubble friction term. Note, in particular, how the Hub- 
ble radius naturally emerges as the scale where the nature of the fluctuating 
modes changes from oscillatory to frozen. 

Considering the equation in a bit more detail, observe that if the equation 
of state of the background is independent of time (which will be the case if 
TL = ip — 0), that then in an expanding background, the dominant mode of 
l|62|) is constant, and the sub-dominant mode decays. If the equation of state 
is not constant, then the dominant mode is not constant in time. Specifically, 
at the end of inflation H < 0, and this leads to a growth of <fi (see the 
following subsection) . 

To study the quantitative implications of the equation of motion (|62f> . 
it is convenient to introduce HO] the variable ( (which, up to correction 
terms of the order V 2 </> which are unimportant for large-scale fluctuations is 
equal to the curvature perturbation 1Z in comoving gauge |41j 1 by 

2(H- 1 6 + 6) , 
C ^+-i— f-JH, 63 
3 1 + w 

where 

w = - (64) 
P 

characterizes the equation of state of matter. In terms of £, the equation of 
motion (|62(l takes on the form 

~(H(l + w) = O + O(V 2 0). (65) 

On large scales, the right hand side of the equation is negligible, which leads 
to the conclusion that large-scale cosmological fluctuations satisfy 

C(l + w) = 0. (66) 
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This implies that except possibly if 1 + w = at some points in time during 
cosmological evolution (which occurs during reheating in inflationary cos- 
mology if the inflaton field undergoes oscillations - see 021 an d [431 l4"4"| for 
discussions of the consequences in single and double field inflationary models, 
respectively) £ is constant. In single matter field models it is indeed possible 
to show that ( = on super-Hubble scales independent of assumptions on 
the equation of state [451 I46| . This "conservation law" makes it easy to relate 
initial fluctuations to final fluctuations in inflationary cosmology, as will be 
illustrated in the following subsection. 



3.5 Application to Inflationary Cosmology 

Let us now return to the space-time sketch of the evolution of fluctuations 
in inflationary cosmology (Figure 1) and use the conservation law 1[66[1 - in 
the form £ = const on large scales - to relate the amplitude of <f> at initial 
Hubble radius crossing during the inflationary phase (at t — ti(k)) with the 
amplitude at final Hubble radius crossing at late times (at t — t/(fc)). Since 
both at early times and at late times <j> — on super-Hubble scales as the 
equation of state is not changing, (|rJrJ)l leads to 

*(*/(*))* ^TO ^W)- (67) 
[1 + w)[ti{k)) 

This equation will allow us to evaluate the amplitude of the cosmological 
perturbations when they re-enter the Hubble radius at time i/(fc), under the 
assumption (discussed in detail in the following section) that the origin of 
the primordial fluctuations is quantum vacuum oscillations. 

The time-time perturbed Einstein equation (the first equation of l|57[)) 
relates the value of at initial Hubble radius crossing to the amplitude of 
the relative energy density fluctuations. This, together with the fact that the 
amplitude of the scalar matter field quantum vacuum fluctuations is of the 
order H, yields 

Cj>{ti{k)) ~ Hy(ti(k)). (68) 

In the late time radiation dominated phase, w = 1/3, whereas during slow-roll 
inflation 

l + w(U(k)) si y(ti(k)). (69) 

Making, in addition, use of the slow roll conditions satisfied during the infla- 
tionary period 

Htpo ~ -V 
H 2 s —£-V, (70) 



we arrive at the final result 
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<Kt f (k)) 



y3/2 



-(U(k)), 



(71) 



V 



which gives the position space amplitude of cosmological fluctuations on a 
scale labelled by the comoving wavenumber fc at the time when the scale 
re-enters the Hubble radius at late times, a result first obtained in the case 
of the Starobinsky model 02] of inflation in 42j, and later in the context of 
scalar field-driven inflation in 15T1 . 

In the case of slow roll inflation, the right hand side of l|71|l is, to a 
first approximation, independent of k, and hence the resulting spectrum of 
fluctuations is scale- invariant. 



4 Quantum Theory of Cosmological Fluctuations 
4.1 Overview 

As already mentioned in the last subsection of the previous section, in many 
models of the very early Universe, in particular in inflationary cosmology, 
but also in the Pre-Big-Bang and in the Ekpyrotic scenarios, primordial 
inhomogeneities emerge from quantum vacuum fluctuations on microscopic 
scales (wavelengths smaller than the Hubble radius) . The wavelength is then 
stretched relative to the Hubble radius, becomes larger than the Hubble 
radius at some time and then propagates on super-Hubble scales until re- 
entering at late cosmological times. In the context of a Universe with a de Sit- 
ter phase, the quantum origin of cosmological fluctuations was first discussed 
in - see for a more general discussion of the quantum origin of fluctua- 
tions in cosmology, and also [513021 for earlier ideas. In particular, Mukhanov 
|48| and Press |53j realized that in an exponentially expanding background, 
the curvature fluctuations would be scale-invariant, and Mukhanov provided 
a quantitative calculation which also yielded the logarithmic deviation from 
exact scale-invariance. 

To understand the role of the Hubble radius, consider the equation of a 
free scalar matter field ip on an unperturbed expanding background: 



The second term on the left hand side of this equation leads to damping 
of ip with a characteristic decay rate given by H. As a consequence, in the 
absence of the spatial gradient term, ip would be of the order of magnitude 
Hip. Thus, comparing the second and the third term on the left hand side, we 
immediately see that the microscopic (spatial gradient) term dominates on 
length scales smaller than the Hubble radius, leading to oscillatory motion, 
whereas this term is negligible on scales larger than the Hubble radius, and 
the evolution of ip is determined primarily by gravity. 



V 2 

ip + 3Hip T if = . 



a 



(72) 
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To understand the generation and evolution of fluctuations in current 
models of the very early Universe, we thus need both Quantum Mechanics 
and General Relativity, i.e. quantum gravity. At first sight, we are thus faced 
with an intractable problem, since the theory of quantum gravity is not yet 
established. We are saved by the fact that today on large cosmological scales 
the fractional amplitude of the fluctuations is smaller than 1. Since gravity 
is a purely attractive force, the fluctuations had to have been - at least in 
the context of an eternally expanding background cosmology - very small in 
the early Universe. Thus, a linearized analysis of the fluctuations (about a 
classical cosmological background) is self-consistent. 

From the classical theory of cosmological perturbations discussed in the 
previous section, we know that the analysis of scalar metric inhomogeneities 
can be reduced - after extracting gauge artifacts - to the study of the evolution 
of a single fluctuating variable. Thus, we conclude that the quantum theory 
of cosmological perturbations must be reducible to the quantum theory of a 
single free scalar field which we will denote by v. Since the background in 
which this scalar field evolves is time-dependent, the mass of v will be time- 
dependent. The time-dependence of the mass will lead to quantum particle 
production over time if we start the evolution in the vacuum state for v. As we 
will see, this quantum particle production corresponds to the development 
and growth of the cosmological fluctuations. Thus, the quantum theory of 
cosmological fluctuations provides a consistent framework to study both the 
generation and the evolution of metric perturbations. The following analysis 
is based on Part II of [H]- 

4.2 Outline of the Analysis 

In order to obtain the action for linearized cosmological perturbations, we ex- 
pand the action to quadratic order in the fluctuating degrees of freedom. The 
linear terms cancel because the background is taken to satisfy the background 
equations of motion. 

We begin with the Einstein-Hilbert action for gravity and the action of a 
scalar matter field (for the more complicated case of general hydrodynamical 
fluctuations the reader is referred to 0) 



where g is the determinant of the metric. 

The simplest way to proceed is to work in a fixed gauge, longitudinal 
gauge, in which the metric and matter take the form 



<p(v>x) = foil) + 8<p(r],x) . 

The next step is to reduce the number of degrees of freedom. First, as 
already mentioned in the previous section, the off-diagonal spatial Einstein 




(73) 



ds 2 = a 2 (T 1 )[(l + 2(t)(r 1 ,x))dT 1 2 - (1 - 2i/>(^ x))dx 2 ] 



(74) 
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equations force ip = <f> since 5T* = for scalar field matter (no anisotropic 
stresses to linear order). The two remaining fluctuating variables </> and (p 
must be linked by the Einstein constraint equations since there cannot be 
matter fluctuations without induced metric fluctuations. 

The two nontrivial tasks of the lengthy jS] computation of the quadratic 
piece of the action is to find out what combination of ip and (j> gives the 
variable v in terms of which the action has canonical form, and what the 
form of the time-dependent mass is. This calculation involves inserting the 
ansatz (|74|l into the action J73J), expanding the result to second order in 
the fluctuating fields, making use of the background and of the constraint 
equations, and dropping total derivative terms from the action. In the context 
of scalar field matter, the quantum theory of cosmological fluctuations was 
developed by Mukhanov |55l Iftfij (see also [S3)- The result is the following 
contribution S^ 2 ' to the action quadratic in the perturbations: 

s w = ly d ^_ v . v . + ^ v 2^ (75) 

where the canonical variable v (the "Mukhanov variable" introduced in [SHj 
- see also [52]) is given by 



v = u I u if T 

with Ji = a! /a, and where 



= a[Sip + ^], (76) 



* = ^. (77) 

In both the cases of power law inflation and slow roll inflation, H. and ip are 
proportional and hence 

z{rj) ~ a( V ) . (78) 

Note that the variable v is related to the curvature perturbation 1Z in comov- 
ing coordinates introduced in 02 and closely related to the variable £ used 
in 39 , 40 : 

v = zK. (79) 
The equation of motion which follows from the action (|75|l is 



v -V 2 v v = , (80) 

z 

or, in momentum space 

// 

vl + k 2 v k - — v k = , (81) 



where v k is the k'th Fourier mode of v. As a consequence of iJJSJ, the mass 
term in the above equation is given by the Hubble scale 
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k 2 H = — ~ H 2 . (82) 



z 



Thus, it immediately follows from Ij81|l that on small length scales, i.e. for 
k > kn, the solutions for Vk are constant amplitude oscillations . These 
oscillations freeze out at Hubble radius crossing, i.e. when k — ku- On longer 
scales (k -C ku), the solutions for Vk increase as z: 

Vk ~ z , k <C fc# ■ (83) 

Given the action (|75|l . the quantization of the cosmological perturbations 
can be performed by canonical quantization (in the same way that a scalar 
matter field on a fixed cosmological background is quantized 

The final step in the quantum theory of cosmological perturbations is to 
specify an initial state. Since in inflationary cosmology, all pre-existing clas- 
sical fluctuations are red-shifted by the accelerated expansion of space, one 
usually assumes (we will return to a criticism of this point when discussing 
the trans-Planckian problem of inflationary cosmology) that the field v starts 
out at the initial time U mode by mode in its vacuum state. Two questions 
immediately emerge: what is the initial time tj, and which of the many pos- 
sible vacuum states should be chosen. It is usually assumed that since the 
fluctuations only oscillate on sub-Hubble scales, that the choice of the initial 
time is not important, as long as it is earlier than the time when scales of 
cosmological interest today cross the Hubble radius during the inflationary 
phase. The state is usually taken to be the Bunch-Davies vacuum (see e.g. 
|58|1. since this state is empty of particles at ij in the coordinate frame deter- 
mined by the FLRW coordinates (see e.g. [SHj for a discussion of this point), 
and since the Bunch-Davies state is a local attractor in the space of initial 
states in an expanding background (see e.g. (HOj). Thus, we choose the initial 
conditions 

(84) 




where here uik = k, and rji is the conformal time corresponding to the physical 
time fj. 

Let us briefly summarize the quantum theory of cosmological perturba- 
tions. In the linearized theory, fluctuations are set up at some initial time ti 
mode by mode in their vacuum state. While the wavelength is smaller than 
the Hubble radius, the state undergoes quantum vacuum fluctuations. The 
accelerated expansion of the background redshifts the length scale beyond 
the Hubble radius. The fluctuations freeze out when the length scale is equal 
to the Hubble radius. On larger scales, the amplitude of Vk increases as the 
scale factor. This corresponds to the squeezing of the quantum state present 
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at Hubble radius crossing (in terms of classical general relativity, it is self- 
gravity which leads to this growth of fluctuations). As discussed e.g. in |rJT] . 
the squeezing of the quantum vacuum state leads to the emergence of the 
classical nature of the fluctuations. 

4.3 Application to Inflationary Cosmology 

In this subsection we will use the quantum theory of cosmological perturba- 
tions developed in this section to calculate the spectrum of curvature fluctu- 
ations in inflationary cosmology. 

We need to compute the power spectrum Vn(k) of the curvature fluctu- 
ation 1Z defined in (|79|l . namely 

K = z~ l v = (j) + Sip— (85) 

The idea in calculating the power spectrum at a late time t is to first relate 
the power spectrum via the growth rate Ij83(l of v on super-Hubble scales to 
the power spectrum at the time tn(k) of Hubble radius crossing, and to then 
use the constancy of the amplitude of v on sub-Hubble scales to relate it to 
the initial conditions l|84|) . Thus 

Vn(k,t) = k 3 TZ 2 k (t) = k 3 z- 2 (t)\v k (t)\ 2 (86) 

= fc3 - 2 ^(l(^|))) 2 '-^«)' 2 
= k 3 z' 2 (t H (k))\v k (t H (k))\ 2 
~ k 3 a- 2 (t H (k))\v k (ti)\ 2 , 

where in the final step we have used (|78|l and the constancy of the amplitude 
of v on sub-Hubble scales. Making use of the condition 

a" 1 (t ff (fc))fc = H (87) 

for Hubble radius crossing, and of the initial conditions l|84|) . we immediately 
see that 

Vn(k,t) ~ k^k^k^H 2 , (88) 

and that thus a scale invariant power spectrum with amplitude proportional 
to H 2 results, in agreement with what was argued on heuristic grounds in 
Section (|3.5|) . 

4.4 Quantum Theory of Gravitational Waves 

The quantization of gravitational waves parallels the quantization of scalar 
metric fluctuations, but is more simple because there are no gauge ambigui- 
ties. Note again that at the level of linear fluctuations, scalar metric fluctua- 
tions and gravitational waves arc independent. Both can be quantized on the 
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same cosmological background determined by background scale factor and 
background matter. However, in contrast to the case of scalar metric fluctu- 
ations, the tensor modes are also present in pure gravity (i.e. in the absence 
of matter). 

Starting point is the action (|73[1 . Into this action we insert the metric 
which corresponds to a classical cosmological background plus tensor metric 
fluctuations: 

ds 2 = a 2 {r/) [dr? - (fly + h i j)dx i dx j ] , (89) 

where the second rank tensor /iy(?7, x) represents the gravitational waves, 
and in turn can be decomposed as 

hij(v, x) = h+(f], x)e+- + h x {rj, x)e* (90) 

into the two polarization states. Here, efj and e|j- are two fixed polarization 
tensors, and h+ and h x are the two coefficient functions. 

To quadratic order in the fluctuating fields, the action separates into sep- 
arate terms involving h+ and h x . Each term is of the form 

iP> = Jd^[h' 2 -(Vh) 2 ], (91) 

leading to the equation of motion 

hi + 2—h' k + k 2 h k = 0. (92) 
a 

The variable in terms of which the action (|91|l has canonical kinetic term is 

Hk = ah k , (93) 

and its equation of motion is 

vl + (k 2 -— )fik = 0. (94) 

This equation is very similar to the corresponding equation (|81|l for scalar 
gravitational inhomogeneities, except that in the mass term the scale factor 
a(r]) is replaced by z(rf), which leads to a very different evolution of scalar 
and tensor modes during the reheating phase in inflationary cosmology during 
which the equation of state of the background matter changes dramatically. 

Based on the above discussion we have the following theory for the gener- 
ation and evolution of gravitational waves in an accelerating Universe (first 
developed by Grishchuk 62 ): waves exit as quantum vacuum fluctuations at 
the initial time on all scales. They oscillate until the length scale crosses the 
Hubble radius. At that point, the oscillations freeze out and the quantum 
state of gravitational waves begins to be squeezed in the sense that 

Vk(v) ~ a(v) > ( 95 ) 

which, from l|93(l corresponds to constant amplitude of h k ■ The squeezing of 
the vacuum state leads to the emergence of classical properties of this state, 
as in the case of scalar metric fluctuations. 
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5 The Trans-Planckian Window 

Whereas the contents of the previous sections are well established, this and 
the following section deal with aspects of the theory of cosmological pertur- 
bations which are currently under investigation and are at the present time 
rather controversial. First, we consider the trans-Planckian issue (this section 
is adapted from [SH]). 

The same background dynamics which yields the causal generation mech- 
anism for cosmological fluctuations, the most spectacular success of infla- 
tionary cosmology, bears in it the nucleus of the "trans-Planckian problem" . 
This can be seen from Fig. 3. If inflation lasts only slightly longer than the 
minimal time it needs to last in order to solve the horizon problem and to 
provide a causal generation mechanism for CMB fluctuations, then the corre- 
sponding physical wavelength of these fluctuations is smaller than the Planck 
length at the beginning of the period of inflation. The theory of cosmologi- 
cal perturbations is based on classical general relativity coupled to a weakly 
coupled scalar field description of matter. Both the theories of gravity and of 
matter will break down on trans-Planckian scales, and this immediately leads 
to the trans-Planckian problem: are the predictions of standard inflationary 
cosmology robust against effects of trans-Planckian physics |64|? 

The simplest way of modeling the possible effects of trans-Planckian 
physics, while keeping the mathematical analysis simple, is to replace the 
linear dispersion relation ui. g = fc p hy S of the usual equation for cosmological 
perturbations by a non standard dispersion relation u). s = ^ phya (k) which 
differs from the standard one only for physical wavenumbers larger than the 
Planck scale. This method was introduced 65, 65] in the context of studying 
the dependence of the thermal spectrum of black hole radiation on trans- 
Planckian physics. In the context of cosmology, it has been shown |67U68ll6*9"| 
that this amounts to replacing k 2 appearing in (|8f \ with k 2 s (n, rj) defined by 



For a fixed comoving mode, this implies that the dispersion relation becomes 
time-dependent. Therefore, the equation of motion of the quantity Vk(r]) takes 
the form (with z(r)) oc a(rj)) 



A more rigorous derivation of this equation, based on a variational principle, 
has been provided [7U] (see also Ref. I7T)). 

The evolution of modes thus must be considered separately in three 
phases, see Fig. 3. In Phase I the wavelength is smaller than the Planck 
scale, and trans-Planckian physics can play an important role. In Phase II, 




(96) 



v 'k + fc cff(M) v k = 0. 



(97) 
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Phase III 



Phase 11 



Phase I 




Fig. 3. Space-time diagram (physical distance vs. time) showing the origin of 
the trans-Planckian problem of inflationary cosmology: at very early times, the 
wavelength is smaller than the Planck scale £p\ (Phase I), at intermediate times 
it is larger than £p\ but smaller than the Hubble radius H^ 1 (Phase II), and at 
late times during inflation it is larger than the Hubble radius (Phase III). The 
line labeled a) is the physical wavelength associated with a fixed comoving scale k. 
The line b) is the Hubble radius or horizon in SBB cosmology. Curve c) shows the 
Hubble radius during inflation. The horizon in inflationary cosmology is shown in 
curve d). 

the wavelength is larger than the Planck scale but smaller than the Hub- 
ble radius. In this phase, trans-Planckian physics will have a negligible effect 
(this statement can be quantified [72j). Hence, by the analysis of the previous 
section, the wave function of fluctuations is oscillating in this phase, 

Vk — B\ exp(—ikri) + B2 cxp^/cry) (98) 

with constant coefficients B\ and Bi- In the standard approach, the initial 
conditions are fixed in this region and the usual choice of the vacuum state 
leads to B\ = l/V^fc, B2 = 0. Phase III starts at the time tu(k) when 
the mode crosses the Hubble radius. During this phase, the wave function is 
squeezed. 

One source of trans-Planckian effects [S3 ||>U on observations is the pos- 
sible non-adiabatic evolution of the wave function during Phase I. If this 
occurs, then it is possible that the wave function of the fluctuation mode is 
not in its vacuum state when it enters Phase II and, as a consequence, the 
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coefficients Bi and B 2 are no longer given by the standard expressions above. 
In this case, the wave function will not be in its vacuum state when it crosses 
the Hubble radius, and the final spectrum will be different. In general, Bi 
and B2 are determined by the matching conditions between Phase I and II. 
If the dynamics is adiabatic throughout (in particular if the a" /a term is 
negligible), the WKB approximation holds and the solution is always given 

by 

Vk ^) ~ /ni /, ^ X PM / fc cffdr) , (99) 
^2fc off (fc,ry) V Jni J 

where rji is some initial time. Therefore, if we start with a positive frequency 
solution only and use this solution, we find that no negative frequency solution 
appears. Deep in Region II where fc c g ~ k the solution becomes 

Wfe ( 7 ?) - ~7kt ex P(-# - ik v), (100) 
V 2k 

i.e. the standard vacuum solution times a phase which will disappear when we 
calculate the modulus. To obtain a modification of the inflationary spectrum, 
it is sufficient to find a dispersion relation such that the WKB approximation 
breaks down in Phase I. 

A concrete class of dispersion relations for which the WKB approximation 
breaks down is 

k 2 cS (k, V ) = k 2 -k 2 \b m \ 

where X(rj) — 2na(r]) /k is the wavelength of a mode. If we follow the evolution 
of the modes in Phases I, II and III, matching the mode functions and their 
derivatives at the junction times, the calculation |^7| 1681 175| demonstrates 
that the final spectral index is modified and that superimposed oscillations 
appear. 

However, the above example suffers from several problems. First, in infla- 
tionary models with a long period of inflationary expansion, the dispersion 
relation H1U1|) leads to complex frequencies at the beginning of inflation for 
scales which are of current interest in cosmology. Furthermore, the initial 
conditions for the Fourier modes of the fluctuation field have to be set in a 
region where the evolution is non-adiabatic and the use of the usual vacuum 
prescription can be questioned. These problems can be avoided in a toy model 
in which we follow the evolution of fluctuations in a bouncing cosmological 
background which is asymptotically flat in the past and in the future. The 
analysis shows that even in this case the final spectrum of fluctuations 
depends on the specific dispersion relation used. 

An example of a dispersion relation which breaks the WKB approximation 
in the trans-Planckian regime but does not lead to the problems mentioned 
in the previous paragraph was investigated in [70]. It is a dispersion relation 
which is linear for both small and large wavenumbers, but has an intermediate 



(101) 
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interval during which the frequency decreases as the wavenumber increases, 
much like what happens in i|l(Jl|) • The violation of the WKB condition occurs 
for wavenumbers near the local minimum of the u(k) curve. 

A justified criticism against the method summarized in the previous anal- 
ysis is that the non-standard dispersion relations used are completely ad hoc, 
without a clear basis in trans-Planckian physics. There has been a lot of recent 
work [23 EH IZ3 HE! on the implication of space-space uncertainty relations 
jZHl ED! on the evolution of fluctuations. The application of the uncertainty 
relations on the fluctuations lead to two effects |^ |H2 ■ Firstly, the equation 
of motion of the fluctuations in modified. Secondly, for fixed comoving length 
scale fc, the uncertainty relation is saturated before a critical time ti(k). Thus, 
in addition to a modification of the evolution, trans-Planckian physics leads 
to a modification of the boundary condition for the fluctuation modes. The 
upshot of this work is that the spectrum of fluctuations is modified. 

In |83|. the implications of the stringy space-time uncertainty relation 



on the spectrum of cosmological fluctuations was studied. Again, application 
of this uncertainty relation to the fluctuations leads to two effects. Firstly, the 
coupling between the background and the fluctuations is nonlocal in time, 
thus leading to a modified dynamical equation of motion (a similar modi- 
fication also results [EE] from quantum deformations, another example of a 
consequence of non-commutative basic physics). Secondly, the uncertainty 
relation is saturated at the time U(k) when the physical wavelength equals 
the string scale l s . Before that time it does not make sense to talk about 
fluctuations on that scale. By continuity, it makes sense to assume that fluc- 
tuations on scale k are created at time ti{k) in the local vacuum state (the 
instantaneous WKB vacuum state). 

Let us for the moment neglect the nonlocal coupling between background 
and fluctuation, and thus consider the usual equation of motion for fluctua- 
tions in an accelerating background cosmology. We distinguish two ranges of 
scales. Ultraviolet modes are generated at late times when the Hubble radius 
is larger than On these scales, the spectrum of fluctuations does not differ 
from what is predicted by the standard theory, since at the time of Hubble 
radius crossing the fluctuation mode will be in its vacuum state. However, 
the evolution of infrared modes which are created when the Hubble radius is 
smaller than l a is different. The fluctuations undergo less squeezing than they 
do in the absence of the uncertainty relation, and hence the final amplitude 
of fluctuations is lower. From the equation i|8fi|) for the power spectrum of 
fluctuations, and making use of the condition 



mug 



Ax phys At > l 2 s 



(102) 



a(U(k)) = kl s 



(103) 



for the time ti(k) when the mode is generated, it follows immediately that 
the power spectrum is scale-invariant 
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V n (k)~k°. (104) 

In the standard scenario of power-law inflation the spectrum is red (Vn(k) ~ 
with n < 1). Taking into account the effects of the nonlocal coupling 
between background and fluctuation mode leads [HSj to a modification of this 
result: the spectrum of fluctuations in a power-law inflationary background 
is in fact blue (n > 1). 

Note that, if we neglect the nonlocal coupling between background and 
fluctuation mode, the result of H 104(1 also holds in a cosmological background 
which is NOT accelerating. Thus, we have a method of obtaining a scale- 
invariant spectrum of fluctuations without inflation. This result has also been 
obtained in |87| , however without a micro-physical basis for the prescription 
for the initial conditions. 

A key problem with the method of modified dispersion relations is the 
issue of back-reaction [%%| 155 ] . If the mode occupation numbers of the fluc- 
tuations at Hubble radius crossing are significant, the danger arises that the 
back-reaction of the fluctuations will in fact prevent inflation. Another con- 
straint arises from the observational limits on the flux of ultra-high-energy 
cosmic rays. Such cosmic rays would be produced |90j in the present Universe 
if Trans-Planckian effects of the type discussed in this section were present. 

An approach to the trans-Planckian issue pioneered by Danielsson [OJ] 
which has recently received a lot of attention is to avoid the issue of the 
unknown trans-Planckian physics and to start the evolution of the fluctuation 
modes at the mode- dependent time when the wavelength equals the limiting 
scale. Obviously, the resulting spectrum will depend sensitively on which state 
is taken to be the initial state. The vacuum state is not unambiguous, and the 
choice of a state minimizing the energy density depends on the space-time 
splitting |22]. The signatures of this prescription are typically oscillations 
superimposed on the usual spectrum. The amplitude of this effect depends 
sensitively on the prescription of the initial state, and for a fixed prescription 
also on the background cosmology. For a discussion of these issues and a list 
of references on this approach the reader is referred to jHS] • 

In summary, due to the exponential red-shifting of wavelengths, present 
cosmological scales originate at wavelengths smaller than the Planck length 
early on during the period of inflation. Thus, Planck physics may well encode 
information in these modes which can now be observed in the spectrum of 
microwave anisotropies. Two examples have been shown to demonstrate the 
existence of this "window of opportunity" to probe trans-Planckian physics 
in cosmological observations. The first method makes use of modified dis- 
persion relations to probe the robustness of the predictions of inflationary 
cosmology, the second applies the stringy space-time uncertainty relation on 
the fluctuation modes. Both methods yield the result that trans-Planckian 
physics may lead to measurable effects in cosmological observables. An im- 
portant issue which must be studied more carefully is the back-reaction of 
the cosmological fluctuations (see e.g. ]P^] for a possible formalism). 
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6 Back-Reaction of Cosmological Fluctuations 

The presence of cosmological fluctuations influences the background cosmol- 
ogy in which the perturbations evolve. This back-reaction arises as a second 
order effect in the cosmological perturbation expansion. The effect is cumu- 
lative in the sense that all fluctuation modes contribute to the change in the 
background geometry, and as a consequence the back-reaction effect can be 
large even if the amplitude of the fluctuation spectrum is small. In this section 
(based on the review we discuss two approaches used to quantify back- 
reaction. In the first approach |961 194) . the effect of the fluctuations on the 
background is expressed in terms of an effective energy-momentum tensor. 
We show that in the context of an inflationary background cosmology, the 
long wavelength contributions to the effective energy-momentum tensor take 
the form of a negative cosmological constant, whose absolute value increases 
as a function of time since the phase space of infrared modes is increasing. 
This then leads to the speculation |971 l9"5] that gravitational back-reaction 
may lead to a dynamical cancellation mechanism for a bare cosmological con- 
stant, and yield a scaling fixed point in the asymptotic future in which the 
remnant cosmological constant satisfies Qa ~ 1< We then discuss (HHl how 
infrared modes effect local observables (as opposed to mathematical back- 
ground quantities) and find that the leading infrared back-reaction contribu- 
tions cancel in single field inflationary models. However, we expect non-trivial 
back-reaction of infrared modes in models with more than one matter field. 

It is well known that gravitational waves propagating in some background 
space-time affect the dynamics of the background. This back-reaction can be 
described in terms of an effective energy- momentum tensor t^ v . In the short 
wave limit, when the typical wavelength of the waves is small compared with 
the curvature of the background space-time, t m „ has the form of a radiative 
fluid with an equation of state p = p/3 (where p and p denote pressure and 
energy density, respectively). As we have seen in previous section, in infla- 
tionary cosmology it is the long wavelength scalar metric fluctuations which 
are more important. Like short wavelength gravitational waves, these cosmo- 
logical fluctuations will contribute to the effective energy-momentum tensor 
r^v . The work of jHEl US] is closely related to work by Woodard and Tsamis 
|l()()l ll()l| who considered the back-reaction of long wavelength gravitational 
waves in pure gravity with a bare cosmological constant. The recent paper 
|9"9"] is related to the work of Abramo and Woodard [1021 1103| who initiated 
the study of back-reaction of infrared modes on local observables. 

We first review the derivation of the effective energy-momentum tensor 
T^y which describes the back-reaction of linear cosmological fluctuations on 
the background cosmology, and summarize the evaluation of this tensor in an 
inflationary cosmological background. This gravitational back-reaction calcu- 
lation is related to the early work on the back-reaction of gravitational waves 
by Brill, Hartle and Isaacson jl()4| . among others. The idea is to expand the 
Einstein equations to second order in the perturbations, to assume that the 
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first order terms satisfy the equations of motion for linearized cosmological 
perturbations discussed in previous section (hence these terms cancel), to 
take the spatial average of the remaining terms, and to regard the resulting 
equations as equations for a new homogeneous metric gfy which includes 
the effect of the perturbations to quadratic order: 



T (o) 

[IIS 



1 fiis 



(105) 



where the effective energy-momentum tensor 7) t „ of gravitational back- 
reaction contains the terms resulting from spatial averaging of the second 
order metric and matter perturbations: 



= < r( 2 ) — <7( 2 ) > 

8irG ' 



(106) 



where pointed brackets stand for spatial averaging, and the superscripts in- 
dicate the order in perturbation theory. 

As analyzed in detail in jHElEl, the back-reaction equation 1)105(1 is co- 
variant under linear space-time coordinate transformations even though t^ v 
is not invariant 2 . In the following, we will work in longitudinal gauge. 

For simplicity, we shall take matter to be described in terms of a single 
scalar field. By expanding the Einstein and matter energy-momentum tensors 
to second order in the metric and matter fluctuations cj> and Sip, respectively, 
it can be shown that the non-vanishing components of the effective back- 
reaction energy-momentum tensor t^ v become 



1 
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8ttG 
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(107) 



and 



a 2 6i 



8ttG 



(24ff 2 + 16F)(0 2 ) + 24H 
4 



+ <(0) 2 )+4<^)--a- 2 ((V0) 2 ) 



;((M 2 



-V"(<po)(5<p 2 ) + 2V'(<po){ct>8<p) 



(108) 



2 See |1()5) . however, for important questions concerning the covariance of the 
analysis under higher order coordinate transformations. 
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where H is the Hubble expansion rate. 

The metric and matter fluctuation variables <p and Sip are linked via the 
Einstein constraint equations, and hence all terms in the above formulas 
for the components of can be expressed in terms of two point functions 
of 4> and its derivatives. The two point functions, in turn, are obtained by 
integrating over all of the Fourier modes of <j>, e.g. 



where <f>k denotes the amplitude of the k'th Fourier mode. The above expres- 
sion is divergent both in the infrared and in the ultraviolet. The ultraviolet 
divergence is the usual divergence of a free quantum field theory and can 
be "cured" by introducing an ultraviolet cutoff k u . In the infrared, we will 
discard all modes k < fcj with wavelength larger than the Hubble radius 
at the beginning of inflation, since these modes are determined by the pre- 
inflationary physics. We take these modes to contribute to the background. 

At any time t we can separate the integral in (jl09[l into the contribution 
of infrared and ultraviolet modes, the separation being defined by setting the 
physical wavelength equal to the Hubble radius. Thus, in an inflationary Uni- 
verse the infrared phase space is continually increasing since comoving modes 
are stretched beyond the Hubble radius, while the ultraviolet phase space is 
either constant (if the ultraviolet cutoff corresponds to a fixed physical wave- 
length), or decreasing (if the ultraviolet cutoff corresponds to fixed comoving 
wavelength). In either case, unless the spectrum of the initial fluctuations 
is extremely blue, two point functions such as (|109fl will at later stages of 
an inflationary Universe be completely dominated by the infrared sector. In 
the following, we will therefore restrict our attention to this sector, i.e. to 
wavelengths larger than the Hubble radius. 

In order to evaluate the two point functions which enter into the ex- 
pressions for t m „, we make use of the known time evolution of the linear 
fluctuations <pk discussed in previous section. On scales larger than the Hub- 
ble radius, and for a time-independent equation of state, (f>k is constant in 
time. From the Einstein constraint equations relating the metric and mat- 
ter fluctuations, and making use of the inflationary slow roll approximation 
conditions we find 




(109) 





(111) 



and 




1 



(112) 



Pbr 



Pbr , 
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The main result which emerges from this analysis is that the equation of 
state of the dominant infrared contribution to the energy-momentum tensor 
Tn V which describes back-reaction takes the form of a negative cosmological 
constant 

Pbr = -pbr w ith Pbr < 0. (113) 

The second crucial result is that the magnitude of pbr increases as a function 
of time. This is due in part to the fact that, in an inflationary Universe, as 
time increases more and more wavelengths become longer than the Hubble 
radius and begin to contribute to pb r . 

How large is the magnitude of back-reaction? The basic point is that since 
the amplitude of each fluctuation mode is small, we need a very large phase 
space of infrared modes in order to induce any interesting effects. In models 
with a very short period of primordial inflation, the back-reaction of long- 
wavelength cosmological fluctuations hence will not be important. However, 
in many single field models of inflation, in particular in those of chaotic 
inflation type inflation lasts so long that the infrared back-reaction 

effects can build up to become important for the cosmological background 
dynamics. 

To give an example, consider chaotic inflation with a potential 

V(<p) = ±m 2 <f?. (114) 

In this case, the values of <fik for long wavelength modes are well known 
(see e.g. 0), and the integral in Hl()9|l can be easily performed, thus yielding 
explicit expressions for the dominant terms in the effective energy-momentum 
tensor. Comparing the resulting back-reaction energy density pb r with the 
background density po, we find 

Pbrjt) 3 m 2 (pl(tj) 

po " 4tt M% 

where Mp denotes the Planck mass. Without back-reaction, inflation would 
end when ipo(t) ~ Mp. Inserting this value into l|115[) . we see that 

if (fio(ti) > ipbr ~ m~ l / 3 Mp , then back-reaction will become important 
before the end of inflation and may shorten the period of inflation. It is 
interesting to compare this value with the scale (pa(U) ~ (p sr — mr 1 ' 
above which the stochastic terms in the scalar field equation of motion arising 
in the context of the stochastic approach to chaotic inflation |ll)7l I1U8| are 
dominant. Notice that since ip sr 3> <fbr (recall that m -C Mp), back-reaction 
effects can be very important in the entire range of field values relevant to 
stochastic inflation. 

Since the back-reaction of cosmological fluctuations in an inflationary 
cosmology acts (see 1)44311 ) like a negative cosmological constant, and since 
the magnitude of the back-reaction effect increases in time, one may speculate 



. <Po[t) 



(445) 
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P7) that back-reaction will lead to a dynamical relaxation of the cosmological 
constant (see Tsamis & Woodard 100 for similar speculations based on the 
back-reaction of long wavelength gravitational waves) . 

The background metric gj^v including back-reaction evolves as if the 
cosmological constant at time t were 

A cS {t) = A + 8nGp br {t) (116) 

and not the bare cosmological constant Aq. Hence one might hope to identify 
<|116|) with a time dependent effective cosmological constant. Since |p& r (t)| 
increases as t grows, the effective cosmological constant will decay. Note that 
even if the initial magnitude of the perturbations is small, eventually (if 
inflation lasts a sufficiently long time) the back-reaction effect will become 
large enough to cancel any bare cosmological constant. 

Furthermore, one might speculate that this dynamical relaxation mecha- 
nism for A will be self-regulating. As long as A e s{t) > 8irGp m (t), where p m (t) 
stands for the energy density in ordinary matter and radiation, the evolution 
of g^u^ is dominated by A c g(t). Hence, the Universe will be undergoing 
accelerated expansion, more scales will be leaving the Hubble radius and the 
magnitude of the back-reaction term will increase. However, once A c tf (t) falls 
below p m {t), the background will start to decelerate, scales will enter the 
Hubble radius, and the number of modes contributing to the back-reaction 
will decrease, thus reducing the strength of back-reaction. Hence, it is likely 
that there will be a scaling solution to the effective equation of motion for 
A c fi(t) of the form 

A cS (t) ~ ^G Pm {t) . (117) 

Such a scaling solution would correspond to a contribution to the relative 
closure density of Qa ~ 1. 

There are important concerns about the above formalism, and even more 
so about the resulting speculations (many of these were first discussed in print 
in |105j ). On a formal level, since our back-reaction effect is of second order in 
cosmological perturbation theory, it is necessary to demonstrate covariance 
of the proposed back-reaction equation l|105f) beyond linear order, and this 
has not been done. Next, it might be argued that by causality super-Hubble 
fluctuations cannot affect local observables. Thirdly, from an observational 
perspective one is not interested in the effect of fluctuations on the back- 
ground metric (since what the background is cannot be determined precisely 
using local observations). Instead, one should compute the back-reaction of 
cosmological fluctuations on observables describing the local Hubble expan- 
sion rate. One might then argue that even if long- wavelength fluctuations have 
an effect on the background metric, they do not influence local observables. 
Finally, it is clear that the speculations in the previous section involve the 
extrapolation of perturbative physics deep into the non-perturbative regime. 

These important issues have now begun to be addressed. Good physical 
arguments can be given |102l 1103) supporting the idea that long- wavelength 
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fluctuations can effect local physics. Consider, for example, a black hole of 
mass M absorbing a particle of mass m. Even after this particle has disap- 
peared beyond the horizon, its gravitational effects (in terms of the increased 
mass of the black hole) remain measurable to an external observer. A similar 
argument can be given in inflationary cosmology: consider an initial localized 
mass fluctuation with a characteristic physical length scale A in an exponen- 
tially expanding background. Even after the length scale of the fluctuation 
redshifts to be larger than the Hubble radius, the gravitational potential as- 
sociated with this fluctuation remains measurable. On a more technical level, 
it has recently been shown that super-Hubblc scale (but sub-horizon-scale) 
metric fluctuations can be parametrically amplified during inflationary re- 
heating |1()9I 1421 1481 14"4"| . This clearly demonstrates a coupling between local 
physics and super-Hubble-scale fluctuations. 

These arguments, however, make it even more important to focus on 
back-reaction effects of cosmological fluctuations on local physical observ- 
ables rather than on the mathematical background metric. In recent work 
|99| . the leading infrared back-reaction effects on a local observable measur- 
ing the Hubble expansion rate were calculated. 

Consider a perfect fluid with velocity four vector u a in an inhomogencous 
cosmological geometry, then the local expansion rate which generalizes the 
Hubble expansion rate H (t) of homogeneous isotropic Friedmann- Robertson- 
Walker cosmology is given by -|<9, where is the four divergence of u a : 

= u%, (118) 

the semicolon indicating the covariant derivative. In the effects of cos- 
mological fluctuations on this variable were computed to second order in 
perturbation theory. To leading order in the infrared expansion, the result is 

= 3^(W+L> 2 )-3^, (HQ) 
a z 2 a 

where the prime denotes the derivative with respect to conformal time. If we 
now calculate the spatial average of 0, the term linear in <p vanishes, and - 
as expected - we are left with a quadratic back-reaction contribution. 

Superficially, it appears from i|119fl that there is a non-vanishing back- 
reaction effect at quadratic order which is not suppressed for super-Hubble 
modes. However, we must be careful and evaluate not at a constant value 
of the background coordinates, but rather at a fixed value of some physical 
observable. For example, if we work out the value of in the case of a matter- 
dominated Universe, and express the result as a function of the proper time 
r given by dr 2 — a(r]) 2 {l + 2<fi)dr] 2 instead of as a function of conformal time 
77, then we find that the leading infrared terms proportional to <p 2 exactly 
cancel, and that thus there is no un-suppressed infrared back-reaction on the 
local measure of the Hubble expansion rate. 

A more relevant example with respect to the discussion in earlier sections 
is a model in which matter is given by a single scalar field. In this case, the 
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leading infrared back-reaction terms in & are again given by (|119fl which looks 
different from the background value ZH . However, once again it is important 
to express in terms of a physical background variable. If we choose the 
value of the matter field ip as this variable, we find after easy manipulations 
that, including only the leading infrared back-reaction terms, 

G{fp) = VS^V^). (120) 

Hence, once again the leading infrared back-reaction contributions vanish, 
as already found in the work of |ll)3j which considered the leading infrared 
back-reaction effects on a local observable different than the one we have 
used, and applied very different methods 3 . 

However, in a model with two matter fields, it is clear that if we e.g. use the 
second matter field as a physical clock, then the leading infrared back-reaction 
terms will not cancel in (9, and that thus in such models infrared back-reaction 
will be physically observable. The situation will be very much analogous to 
what happens in the case of parametric resonance of gravitational fluctuations 
during inflationary reheating. This process is a gauge artifact in single field 
models of inflation 021 ( see a ls° |45Ulll)l|4l)| h but it is real and unsuppressed 
in certain two field models 031^5- I* 1 the case of two field models, work on 
the analysis of the back-reaction effects of infrared modes on the observable 
representing the local Hubble expansion rate is in progress. 

Provided that it can indeed be shown that infrared modes have a nontriv- 
ial gravitational back-reaction effect in interesting models at second order in 
perturbation theory, it then becomes important to extend the analysis beyond 
perturbation theory. For initial attempts in this direction see |1 1 II 1112] . 
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